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The variances of the fluid-particle acceleration and of the pressure-gradient and viscous force 
are given. The scaling parameters for these variances are velocity statistics measureable with a 
single-wire anemometer. For both high and low Reynolds numbers, asymptotic scaling formulas are 
given; these agree quantitatively with DNS data. Thus, the scaling can be presumed known for all 
Reynolds numbers. Fluid-particle acceleration variance does not obey K41 scaling at any Reynolds 
, number; this is consistent with recent experimental data. The non-dimensional pressure-gradient 

O ■ variance named \t/Xp is shown to be obsolete. 



I. INTRODUCTION 



Accelerations in turbulent flow are violent and are important to many types of studies (La Porta et al, 2001). 
The accelerations discussed here are caused by viscosity, v\j\u. L , the pressure gradient, —d Xi p, and the fluid-particle 



' acceleration, a-i, which are related by the Navier-Stokes equation: 

^ . a-i = Dm/Dt = -d Xi p + vd Xn d Xn Ui, 

where Dui/Dt denotes the time derivative following the motion of the fluid particle, and Ui is velocity. Here, v is the 
, kinematic viscosity, = d Xn d Xn is the Laplacian operator, d denotes differentiation with respect to the subscript 
variable, summation is implied by repeated indices, p is pressure divided by fluid density; density is constant, e is 
energy dissipation rate per unit mass of fluid. 

Kolmogorov's (1941) scaling (K41 scaling) uses e and v as parameters, and is based on the assumption that 
local isotropy is accurate for a given distance r between two points of measurements, x and x' = x + r, r = |r|. 
I— — '. Refinement of K41 to include the effects of turbulence intermittency (Kolmogorov, 1962) leads to quantification of 
the Reynolds number dependence of the deviation from K41 scaling. The choice of Reynolds number has been 

1/2 

R\ = u rms Ar/V, based on Taylor's (1935) scale At = u rms / ^(<9 Xl iti) 2 ^) and u rms = (ufj 1 ^ 2 , where subscript 1 

denotes components along the axis parallel to r; also, from e = I5v ^(9 Xl ui) 2 ^, R\ — u 2 rms / (eis/15) 1 ^ 2 . Quantities 

, at x' are denoted u[, p' , and A«, = Ui — u[, Ap = p — p' ', etc.. 
t— ( ■ During 1948-1951 [Heisenberg (1948), Obukhov & Yaglom (1951), Batchelor (1951)], the pressure-gradient 

correlation, ( d Xi p (x, t) d x >.p (x', t)\ mean-squared pressure gradient, (d Xi pd Xi p), and, closely related to those, the 



pressure structure function, Dp (r) = (^(Ap) 2 J were related to the velocity structure function D\\ (7-) = ^(Aui) 2 

,£h by means of the assumption that Ui and u[ are joint Gaussian random fields. In fact, the essential approximation, 
now known to be poor, is that Aui must be Gaussian (Hill, 1994). One result of that joint Gaussian theory is that 
j^H *' pressure-gradient acceleration has as unobservably small effect from intermittency as does D\\. In fact, pressure- 
gradient statistics are strongly affected by intermittency [Hill & Wilczak (1995, hereafter HW), Hill & Thoroddsen 
(1997, hereafter HT), Hill & Boratav (1997, hereafter HB), Vedula & Yeung (1999, hereafter VY), Gotoh & Rogallo 
(1999, hereafter GR), Gotoh & Fukayama (2001, hereafter GF), Nelkin & Chen (1998), Antonia et al (1999), La 
Porta et al. (2001)]. 



A advanced theory (HW) relates Dp(r), (d Xi p (x, t) d x >.p (x', t)J , (d Xi pd Xi p), and the pressure spectrum to 

the fourth-order velocity structure function: DyH (r) = (AuiAuj AukAui). The theory allows calculation of such 
pressure-gradient-related statistics from components of T>ijki and therefore by means of hot-wire anemometry, as in 
HT. This theory is valid for all Reynolds numbers and is based on local isotropy without further assumptions. The 
degree to which the theory's predictions are accurate must depend on how anisotropic the large scales are and how 
large the Reynolds number is; isotropy must be approached as Reynolds number becomes small. A advanced theory for 
the correlation (v\/'^Uiv\7'l i ,v! j ) is given by HT; HT's relationship of this correlation to the third-order velocity struc- 
ture function Dyj. (r) = (AmAujAuk) has advantages (HT) over the 1948-1951 theory that related (v\7^.Uii'\7'^,u' i ) 
to Dn. The advanced theories have been used to compare Dp (r) calculated from T>ijki (r) with Dp (r) calculated 

from DNS pressure fields (figure 1 of HB), as well as the corresponding calculations for (d Xi p (x, t) d x /.p (x', t) j and 
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(v\/xUiL'\/xiu' i ) (figures 12, 13 of VY). Since the advanced theories use only the Navier-Stokes equation, incompress- 
ibility, and local isotropy, comparisons of data with the theory give a measure of the local anisotropy of the data, of 
numerical limitations, or of inaccuracy of Taylor's hypothesis (when used). 

Most studies of turbulent acceleration use the traditional approach of determining the R\ dependence that 
results from use of K41 scaling of acceleration statistics [e.g., VY, GR, GF, Antonia et al. (1999), La Porta et al. 
(2001)]. The resultant deviation from K41 scaling, i.e. the R\ dependence, is often called 'anomalous'. There is 
no anomaly when the advanced theory is employed. Because R\ contains u rmsi it is affected by the large scales 
where anisotropy is possible, and R\ is therefore not a parameter of the advanced theory. However, to compare the 
advanced theory with the existing body of empirical knowledge, the advanced-theory's scales must be expressed in 
terms of K41 scales, thereby producing dependence on R x and e. However, K41 scaling parameters and R\ are not 
the scales within the advanced theory. 



II. SCALING OF MEAN-SQUARED PRESSURE GRADIENT 

For locally isotropic turbulence, HW gave the relationship between the mean-squared pressure gradient and 
the fourth-order velocity structure function: 

oo 

(d Xi pd XiP ) = X = 4 J r- 3 [D lln (r) + D aaaa (r) - 6D 110(3 (r)] dr, (1) 
o 

where x is shorthand for the mean-squared pressure gradient for the case of local isotropy. In ([!]) , -Dim (f) , D aaaa (r) , 
and Dupp (r) are components of T)ijki (r); a and (3 denote the Cartesian axes perpendicular to r, and the 1-axis is 
parallel to r. Repeated Greek indices do not imply summation. The result ([!]) applies for all Reynolds numbers and 
without approximation other than local isotropy. 

Defining H x as the ratio of the integral in (Q) to its first term, HW wrote (0) as 

X = 4tf x / r- 3 D uu (r) dr. (2) 



o 



Equivalently, @ defines H x . The purpose of (|J) as stated in HW is that if the Reynolds number variation of 
H x is known, then (Q) enables evaluation of x by calculating the integral in (^) using data from a single hot-wire 
anemometer. Further, HW argued that H x is a constant at large Reynolds numbers. VY evaluated H x by means 
of DNS data and found that it is constant at a value of about 0.65 for 80 < R\ < 230, 230 being their maximum R\. 
Their H x only decreased to about 0.55 at R x — 20. Hill (1994) gives H x — ► 0.36 as R\ — > [on the basis that the 
joint Gaussian assumption can be used in this limit and by use of a formula for the velocity correlation for R\ — ► 
given by Batchelor (1956)]. 

It is useful to express the integral in (Q) in terms of quantities that have been measured in the past. On the 
basis of empirical data described in Appendix A, the approximation for high Reynolds numbers is 

X ^ 3.1H x e 3/2 iy- 1/2 F - 79 ~ 3.9H x e 3/2 iy- 1/2 R° x 25 for i? A > 400, (3) 

where F = (^(d Xl ui) 4 ^ / (^(d Xl ui) 2 ^ is a velocity-derivative flatness. For H x — 0.65, (^) agrees quantitatively with 

the DNS data in Table 1 of GF for R\ > 387; thereby, the estimated limitation supported in Appendix A, i.e., 
ft a ^ 10 3 , seems too conservative. That is why the limitation R\ > 400 is given in (||). 

The case of low Reynolds numbers is in Appendix B where Taylor's scaling and data from VY are used; the 
result is 



X 



0.11e 3/2 is- 1/2 R x for R x < 20, (4) 



which is shown in figure 1 of VY. Neither (|]) nor (|3|) is K41 scaling because of their R\ dependence. 

The data of Pearson & Antonia (2001) reveal how the approximation (^) is approached as R\ increases. The 
inner scale of -Dim (?*) is denoted by £ and is defined in Appendix A as the intersection of viscous- and inertial-range 
asymptotic formulas for £>im (r). Thus, £ is a length scale in the dissipation range. In figures 4 and 5 of Pearson 
& Antonia (2001), the scaled components of Dim (r) increase most rapidly at r > £ as R\ is increased until an 
inertial range is attained. This implies that the integral in (||) will approach the asymptote (||) from below. That 



2 



has been observed, as shown in figure 1, wherein the DNS data of VY and GF are plotted with the asymptotic 
formulas (|) and (Q); those asymptotes are graphed to R\ = 100 and 40, respectively. Because H x = 0.65 was 
used in (|) to obtain figure 1, and because VY found H x ~ 0.65 for 80 < R\ < 230, it appears that H x remains 
constant at about 0.65 for R\ > 80. There was no adjustment of (|J) to cause agreement with the DNS. The 
agreement is surprising because the empirical data used in Appendix A suggests at least a 15% uncertainty of the 
coefficients in (||). As a reminder of this fortunate circumstance, a ±15% error bar is shown at R\ = 10 3 in figure 

1. Previously, (d Xi pd Xi p) /e 3 / 2 ^ -1 / 2 oc R 1 ^ 2 has been reported (VY, GR); this is a good fit to the data of VY in the 
range 130 < R\ < 300, and (d Xi pd Xi p) /e 3/2 ^ 1/2 cc R° x 62 is a good fit to the data of GR in the range 39 < R\ < 170. 
SEE THE FIRST PAGE IN ANCILLARY POSTSCRIPT FILE 

FIG. 1. Figure 1. K41 scaled mean-squared pressure gradient (sum of squares of the 3 components). VY data: asterisks; 
GF data: triangles; equations (2.3, 2.4): lines. 15% error: vertical bar. 

In comparison with figure 1, those power laws are local fits over limited ranges of Rx', they are not asymptotic power 
laws. On the basis of their equation (59b), HW gave the first prediction of the increase of (d Xi pd Xi p) /e 3 / 2 ^ -1 / 2 with 
Reynolds number; (H) is a refinement of equation (59b) of HW. GR noted that equation (59b) gave a weaker 

1/2 

dependence on R\ than their observed ~ RJ dependence, found from DNS for which R\ < 175. The reason is now 
apparent from figure 1: R\ < 175 is too low to use (g) or equation (59b) of HW. Figure 1 does not support the 
multifractal result (i.e., i?° 135 , about half the slope of (2.3)] given by Borgas (1993). 



III. SCALING OF v 2 {pljJi 

Derivation of v 2 ^|V 2 u;| 2 ^) from the spatial correlation of ^V 2 it; is given by Hill (2001). It suffices here to 
state the various equivalent formulas: 

Vii (0) = v 2 (|V 2 Ul | 2 ) = 12, 2 ((d 2 Xl u,) 2 ) = 35, 2 ((d 2 XlUl ) 2 ) , (5) 

where d 2 l = d 2 /dx 2 , and for local stationarity, 

4 

Vu (0) = v {ujiWjSij) = --v (stjSjkSki) (6) 
= -™v ((d Xl u p f d XlUl ) = -^(( dxiUl f) = O.SOe 3 / 2 ^ 1 / 2 \S\ , (7) 

x3\ ,// Q ^2\ 3 / 2 



where a;, is vorticity and s^- is the rate of strain, and S = ((c? xi iti) y / (^(d Xl ui) J is the velocity-derivative 

skewness. |5| is known (Sreenivasan & Antonia, 1997) to increase with increasing R\. Thus, all statistics in (H) and 
(^) have the same increase with increasing R\ when they are nondimensionalized using K41 scaling. On the other 
hand, S is approximately constant over the range of about 20 < R\ < 400 (Sreenivasan & Antonia, 1997) such that 
Vu (0) approximately follows K41 scaling in that range. 

The most carefully selected data for S and F at high Reynolds numbers are those of Antonia et al. (1981), 
which data are in agreement with data at R\ = 10 4 by Kolmyansky, Tsinober & Yorish (2001). The data of Antonia 
et al. (1981) are used for F in (g) and give \S\ ~ 0.5 (i? A /400) 11 for R x > 400. Substituting |5| ~ 0.5 (i? A /400) 11 
in gives 



Vu (0) ^ 0.08i? A for Rx > 400. (8) 

For R x < 20, Tavoularis, Bennett & Corrsin (1978), Herring & Kerr (1982) and Kerr (1985) show that \S\ decreases 
and does so more rapidly as Rx — * 0. The data of Herring and Kerr (1982) suggest that \S\ ~ R\/5 for Rx < 1. 
Although stronger empirical evidence would be helpful, |<S| ~ Rx/5 will serve as the asymptotic formula for Rx < 1, 
in which case (|7j) becomes 

Vu (0) ^ 0.06e 3/2 ^" 1/2 i? A for R x < 1. (9) 
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IV. SCALING OF THE MEAN-SQUARED FLUID-PARTICLE ACCELERATION 



For locally isotropic turbulence, the mean-squared fluid-particle acceleration An (0) is x + Vu (0) because the 
correlation of V 2 ,?4 anc » d x p vanishes by local isotropy (Obukhov & Yaglom, 1951). Thus, for any Reynolds number 
for which local isotropy is valid 

oc 

A u (0) = AH X J r- 3 D nil (r) dr - ^ v ((d XlUl f) . (10) 
o 

Therefore, An (0) scales with the sum of two terms that behave differently with Reynolds number. 
Use of (|) and (^ |) in (0) gives for R x > 400, 

An (0) a e^ 2 v- 1 ' 2 (2.QF 079 + 0.3 |5|) ~ e 3 / 2 ^ 1 / 2 (2.5R° X 25 + 0.08i?°- n ) . (11) 

At i?A = 400, the term from \ is 70 times greater than the term from Vu (0), and x increases much faster than Vu (0) 
for increasing R\. The theory of 1948-1951 used the joint Gaussian assumption and thereby greatly underestimated 
(dxiPdxiP) (HW, HT). That theory gives equation (3.18) of Obukhov & Yaglom (1951), which is the same as ( p"l| ) 
with the exception that 1.1 appears in place of 2. OF 79 . Not only is the magnitude of 1.1 {S^ 1 smaller than 

2. OF 79 by a factor of 5 at R\ — 400, in addition, |5| decreases with further increases of R\ contrary to the 
increase of 2.0F - 79 . An empirical result that seems accurate for a variety of flows for R\ > 400 (Champagne, 1978; 
Antonia et at, 1981) is |5| = 0.25F 3 / 8 , such that Q) can be written as A u (0) ~ e 3 / 2 ^ 1 / 2 (2.0F 079 + 0.075F - 375 )' 
for R x > 400. 

For low Reynolds numbers, use of (0) in ( ^Tj| ) gives 

^(0)~£ 3 / 2 j/- 1 / 2 (0.11i? A + 0.3|S'|) fori? A <20. (12) 

From figure 8 of Herring & Kerr (1982) one sees that 0.11J?a > 0.3 |5| even at their minimum R\ of 0.5, and that 
0.11i?A increases rapidly relative to 0.3 \S\ as R\ increases. Thus, I fi2\ ) shows that the term from x 1S the larger 
contribution to An (0) for all R\. The behaviour of An (0), x, and Vu (0) for moderate R\ is shown particularly well 
in figure 1 of VY. 



V. DISCUSSION 



A. Obsolescence of Xt/Xp 

I 2V 1 / 2 / 2V 1 / 2 

The length scales Ap = u 2 ms / ((d Xx p) j and At = u rms / ((d Xl ui) ) were introduced by Taylor (1935) 

(he included a factor of \[2 that has historically been dropped from these definitions), and the ratio At/Ap, and he gave 

the first evaluation of At/Ap from turbulent diffusion measurements. (At/Ap) 2 = {[d Xl p) 2 ^ / u 2 ms (^(dx^i) 2 ^ 

is a scaled mean-squared pressure gradient; that scaling depends on a large-scale parameter, u rrns . As such At/Ap 
is not relevant in the advanced theory except in the limit of R\ — * 0. Batchclor (1951) obtained, for very large 

Reynolds numbers a dependence of Xt/Xp oc R x , which is not correct because it is based on the joint Gaussian 
assumption. By attempting evaluation of At/Ap, Batchelor (1951) was, in effect, attempting to enable determination 
of x from measurements of velocity variance and energy dissipation rate. The advanced theory replaces (At/Ap) 
with H x . Evaluation of H x allows x to be determined from measurement of a single velocity component and the 
simple formula (||). Figure 3 of GR shows At/Ap versus R\ and reveals the following: i) the strong R\ dependence 
of At/Ap that H x does not have; ii) the deviation of the R\ dependence of Xt/Xp from that predicted by the joint 
Gaussian assumption; (this was also found by VY) and iii) R\ < 20 is required to approach the low-Reynolds-number 
asymptote. Whereas H x depends only on the small scales of turbulence, the dependence of Xt/Xp on the large scales 
via u rms shows that Xt/Xp is not relevant in the advanced theory of local isotropy, except in the limit of R\ — > 0; 
for that limit HW shows that H x oc (Xt/Xp) 2 . 
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B. Acceleration data 



Pioneering technology for measuring turbulence accelerations is being developed at Cornell. La Porta et al. 
(2001) report fluid-particle acceleration measured in a cylindrical enclosure containing turbulent water driven by 
counter-rotating blades. Here, 'x' and l y' axes are transverse and parallel to their cylinder axis, respectively. The 
acceleration's flatness factor (a*) / (a 2 ,) in their figure 3 reaches a maximum at R\ « 700, and is decreased at their 
next-higher R\ value, namely 970, the same is true for their K41-scaled acceleration variances (a 2 ) /e 3 / 2 ^ -1 / 2 in their 
figure 4, where i = x and y ). From (|TT|), An (0) /e 3 / 2 ;/ -1 / 2 = 3 (a 2 ) /e 3 / 2 ;/ -1 / 2 is monotonic with R\, unlike in 
figure 4 of La Porta et al. (2001). This suggests that the cause of the maxima in flatness and variance have the same 
cause. Their estimates of R\ and s both depend on the choice of a velocity component; because the turbulence is 
anisotropic, the choice of another velocity component will shift their data points along both ordinate and abscissa. 
This nonuniversality of their scaling is implicated by the disappearance of the maximum in the flatness of d Xl ui as 
presented by Belin et al. (1997) when R\ is replaced by a universal Reynolds number (Hill, 2001). Belin et al. 
(1997J measure near one counter-rotating blade whereas LaPorta et al. (2001) measure in the flow's center. The two 
cylinders have different aspect ratios. It is nevertheless instructive to substitute the values of F and IS*! measured by 
Belin et al. (1997) into (pi]), divide by 3 to obtain the variance of one component of acceleration, and compare the 
result with the data of La Porta et al. (2001 ). This is done in figure 2, wherein the data of VY and GF are shown 
to agree with the data of Belin et al. (1997). For the Belin et al. data, ( pi] ) produces a maximum similar to that of 
the La Porta et al. (2001 ) data. 

SEE THE SECOND PAGE IN ANCILLARY POSTSCRIPT FILE 

FIG. 2. Figure 2 K41 scaled mean-squared acceleration component. La Porta et al. (2001) data, x-component: asterisks; 
y-component: squares; Belin et al. (1997) data in equation (pi]): diamonds; combined DNS data of VY and GF: triangles; 
eq.(4.3) divided by 3: solid line. 

As Belin et al. (1997,) point out, the maximum in their data for F might be specific to the flow between counter- 
rotating blades. If so, the same is likely true of the data of La Porta et al. (2001) such that their data do not support 
K41 scaling of acceleration, and therefore do not contradict the advanced theory. Another possibility is that the data 
of LaPorta et al. (2001) at R\ = 970 is underestimated for unknown reasons. The conclusion suggested by figure 
2 and the above uncertainties in interpretation of the data is that the data supports the scaling given here and that 
such important acceleration measurements must continue. 



VI. CONCLUSION 

The asymptotes (0) and (|]) combined with the DNS data in figure 1 determine (d Xi pd Xi p) for all Reynolds 
numbers. For all Reynolds numbers the advanced-theory scaling is that (d Xi pd Xi p) scales with the integral in ([!]). For 
R\ < 20, (d Xi pd Xi p) scales with s 3 / 2 v~ 1 / 2 R\. Because H x is apparently constant for R\ > 80, (d Xi pd Xi p) scales with 

oo 

J r~ 3 L>im (r) dr [i.e., the integral in (Q)], and (d Xi pd Xi p) scales approximately with e 3 / 2 ^ _1 / 2 F 1 °- 79 for R\ > 400, 
o 

and e 3 / 2 ^ _1 / 2 F - 79 is a good approximation when R\ as small as 200 (see figure 1). Given that H x is constant 
for R\ > 80, (d Xi pd Xi p) could be obtained for R\ > 80 using data from a single-wire anemometer by evaluating the 
integral in (Q); DNS is not necessary. Using velocity data, it is more accurate to evaluate the integral in (|^) than 
its approximation e 3 / 2 j/ _1 / 2 i^ - 79 because evaluating F requires greater spatial resolution than does evaluation of the 
integral for the same level of accuracy. Evaluating H x from (^) using DNS at R\ > 230 would be useful. 

Now, Vn (0), does scale with any of the derivative moments in (||). It does not scale as in the K41 prediction 
(i.e., £ 3 / 2 z; -1 / 2 ) except for those R\ at which S is constant. The statement in VY that Vn (0) does obey K41 scaling 
is based on their data, which are within the R\ range where S is constant. 

Fluid-particle acceleration variance, An (0), does not scale as in the K41 prediction (i.e., e 3 / 2 ^ -1 / 2 ) at large 
Reynolds numbers because of the factor (2.5i?." 25 + 0.08i?° n ) in (pl|). An (0) does not approach K41 scaling as 
R\ — > because |]) and (Q give An (0) ~ 0.17e 3/2 ^~ 1/2 i?A for R x < 1. For all Reynolds numbers, fluid particle 
acceleration does scale with the sum of velocity statistics that appears on the right-hand side of ([To]) . 

The advanced theory is devoid of statistics of the large scales. It seems paradoxical that R\, which depends 
on the large scales through u rmsi is used above to delineate asymptotic regimes. However, an alternative Reynolds 
number that depends only on small scales (Hill, 2001) makes the advanced theory self-contained. Use of existing 
phenomenology caused both e and R\ to appear in this paper. However, practical applications result. Turbulent 
acceleration-induced coalescence of droplets might be key to understanding rain initiation from liquid-water clouds 
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(Shaw & Oncley, 2001). Radars can measure u rms and e, then R\ = u 2 ms / (ev/15) 1 ^ 2 can be determined; then the 
three acceleration variances can be determined from equations given here. The present results thereby support radar 
remote sensing of clouds and cloud microphysical research. 

Appendix A: High-Reynolds-number asymptote 

The lognormal model of Kolmogorov (1962) is used here; the result is found to be insensitive to the in- 
termittency model used. The inertial range formulas are: Dim ( r ) = C"£ 4 ' 3 r 9 L 2 ' 1 ' 9 , q = (4/3) — (2/i/9), and 
D\i (r) = Ce 2 / 3 r p L~''/ 9 , p = (2/3) + (/i/9); L is the integral scale; [i = 0.25 is used (Sreenivasan & Kailasnath, 

1993), as is C — 2 (Sreenivasan, 1995). Viscous-range formulas are used; they are: %! (r) = ^(<9 xi ui) 4 ^) r 4 

andL>n(r) = ^(^ui) 2 ^ 2 = (e/15^)r 2 . The inner scale of -Dim(r), named is defined by equating the 
inertial-range formula with the viscous-range formulas at r = I. In the integrand in (Q), D\m i r ) can be 
scaled by Dim (!) and r by t. Doing so, the integral equals (3/2) I 2 (^(d Xl ui) 4 ^ [HW showed that the re- 
maining dimcnsionless integral has a value of 3/2 for large Reynolds numbers; this is based on use of an equa- 
tion for Dim (f) that is the same as equation (12) by Stolovitzky, Sreenivasan & Juneja (1993), who demon- 
strate its empirical basis]. Next, the otherwise irrelevant Taylor's scale At is introduced to make use of pub- 
lished empirical data. The definition of i and the inertial range formulas are used to obtain {£/rf) A q = 

(15C) 2 (A!M) 2M/9 F(A T )/F; F(At) = D nn (At) / [Dn (X T f; F = ((d Xl ui) 4 ) / ((^ lUl ) 2 ) 2 . The essen- 
tial approximation is that At is in the inertial range (Antonia et al. 1982; Pearson & Antonia, 2001). Then, 
(3/2) ^(c^r) 4 ) = (3/2)C 4 /( 4 -«)l5 2 («- 2 >/( 4 -«) [(X 2 T /r,L) 2tl/9 F (\ T ) / Fe 3 / 2 ^ 1 / 2 . For 10 2 < R x < 

5 x 10 3 , Zocchi et al. (1994) have \ 2 T /Lr) = 30i?^ 1/2 such that (A 2 ,/Lt;) (2ai/9)(2/(4_9)) = (i? A /900)" 2Al/[9(4 - 9)1 ; 
the exponent of X\jhr\ is about 0.04; so R\ can vary greatly but (Af^/L^) 00 remains near unity. The most con- 
sistent data for F at high Reynolds numbers are those of Antonia et al. (1981), which give F ~ 1.36i2?' 31 ; the 
same data (Antonia et al., 1982) give F (At) cc and F (At) ^ 7.1 at R x = 9400 so F (At) ^ 0.72i?^. Thus, 

(3/2) t 2 {{d Xl uif) ~ 0.714^ - 064 Fe 3 / 2 ! /- 1 /2 = 0.76F°- 79 e 3 / 2 iy~^ 2 = 0.97< 25 e 3 / 2 ^- 1 / 2 , i.e., 



r- 3 Dmi (r) dr ~ Q.nF 01 ^ 3 ' 2 ^ 1 ' 2 ~ Q.^R^e 3 ' 2 ^ 1 ' 2 . (Al) 



Given that [i ~ 0.25, (Al) is insensitive of the value of fj, and is therefore also insensitive to the choice of intermittency 



model. Shaw & Oncley (2001) used data from the atmospheric surface layer at R\ — 1500 to obtain that (Al) 
balances to within the accuracy of their value of F, i.e., about 15%. 

The data of Pearson & Antonia (2001) show Dpppp (At) /Dim (At) becoming constant as R\ — > 10 3 from 

oo 

below; such a constant value is a reasonable criterion for the integral in (|l|) to be proportional to J r~ 3 Dim (r) dr as 

o 

in (|2|). They show the variation of Dim (r) and Dg ggg (r) as R\ varies from 38 to 1200 such that an inertial range 



appears at the larger R\. The above approximation (Al) requires that the integral in (Al) approximately converges 
at its upper limit within the inertial range, for which a reasonable criterion is that there be about one decade of the 
power-law. The data of Pearson & Antonia (2001) show such an extent of the power law as R\ ~ 10 3 is attained. 
Thus, R\ ~ 10 3 is a well-supported lower bound for the high-Reynolds-number approximation (|Al|) . 
Appendix B: Low- Reynolds-number asymptote 

The joint Gaussian assumption is not used here. Taylor's (Taylor, 1935) scaling is used; i.e., scales At and 
u rms are used when R\ — > 0. Taylor's scaling gives Dim {f) oc w 4 ms , H x approaches a constant as R\ — > 0. 

oc 

Let x = t/Xt- Then (Q) can be written as x « u tms^T 2 I x ~ 3 Dmi(x) /uf ms dx. The dimensionless integral 

o 

is a number, so x a w 4 ms Ay 2 ; hence xl {s 3 l 2 v~ x l 2 ^ oc R\. This behaviour is shown in figure 1 of VY where it 
appears to become accurate between R\ = 20 and 40. In their Table II, xl (3e 3 / 2 ^ -1 / 2 ) = 0.74 at R\ = 21. 
Thus, x — 0.106e 3 / 2 ^ _1 / 2 i?A for R\ < 20. Compare 0.106 with the prediction of the joint Gaussian assumption: 
6/15 3 / 2 = 0.103 (Hill, 1994). 
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